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Atomic carbon in PSS 2322+1944, a quasar at redshift 4.12 

J. A. BeeW, P. Cox^, D. Downes\ A. Omont^, F. Bertoldi^ C.L. Carilli^ 



o 
o 

(N 
> 

o 



> 
in 

00 

o 



o 



IRAM, 300 rue de la Piscine, 38406 St-Martin-d'Heres, France 
LERMA, Observatoire de Paris, F-75014 Paris, France 

Institut d'Astrophysique Spatiale, Universite de Paris Sud, F-91405 Orsay, France 

Institut d'Astrophysique de Paris, CNRS and Universite de Paris VI, 98b bd. Arago, F-75014 Paris, France 
Max-Planck-Institut fiir Radioastronomie, Auf dem Hugel 69, D-53121 Bonn, Germany 
National Radio Astronomy Observatory, P.O. Box, Socorro, NM 87801, USA 



Received September 21, 2004 / Accepted October 27, 2004 

Abstract. We report the detection of the fine-structure line of neutral carbon in the z = 4.12 quasar 

PSS 2322-1-1944, obtained at the IRAM Plateau de Bure interferometer. The [Cl] ^Pi Po line is detected 
with a signal-to-noise ratio of ~ 6 with a peak intensity of « 2.5 mjy and a velocity-integrated line flux of 
0.81 ± 0.12 Jykms"^. Assuming an excitation temperature of 43 K (equal to the dust temperature), we derive a 
mass of neutral carbon (corrected for magnification) of Mci ~ 1.2 x 10^ M©. In PSS 2322-1-1944, the cooling due 
to C is about 6 times smaller than for CO, whereas the CO and C cooling represents ~ 10~* of the far- infrared 
continuum and more than half of the cooling due to C*". 

Key words, galaxies: formation - galaxies: starburst - galaxies: high-redshift - quasars: emission lines - quasars: 
individual: PSS 2322-1-1944 - cosmology: observations 
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1. Introduction 

With the detection of dust and molecular gas in sources at 
high redshift, it has become possible to probe the physical 
conditions of the interstellar medium in galaxies and in 
the hosts of quasi-stellar objects (QSOs) at cosmological 
distances. The high-z sources detected in ^^CO (to date 
30 sources between 1.44 < z < 6.42 - see, e.g., Greve et al. 
2004) have massive reservoirs of warm and dense molecu- 
lar gas (a few 10^" — lO^'^ Mq), which are predominantly 
excited by extreme starbursts with implied star formation 
rates « lO^MQyr"^. Multiline CO studies are available 
in only a few cases, and the detection of species other than 
carbon monoxide is reported in only a couple of sources. 
A remarkable example is the Cloverleaf, a strong gravita- 
tionally lensed QSO at z = 2.56 where four CO transitions 
were detected, together with the two fine-structure lines 
of neutral carbon [Ci] (Barvainis et al. 1997; Weifi et al. 
2003) and the J=l->0 transition of HON (Solomon et al. 
2003). 

Atomic carbon is an important probe of the neutral 
dense gas. It is a good tracer of molecular gas in external 
galaxies and plays a central role in the cooling of the gas 
(Gerin & Phillips 1998, 2000). In high-z sources, the detec- 
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* Based on observations obtained with the IRAM Plateau de 
Bure interferometer. 



tion of the [Cl] lines enables to obtain further constraints 
on the physical conditions of the interstellar gas in addi- 
tion to those obtained from the CO transitions (Weifi et 
al. 2003), providing useful information on the gas column 
density, the thermal balance, and the UV illumination. 

PSS 2322-1-1944 is an optically luminous, gravitation- 
ally lensed QSO a,t z — 4.12 which was studied in detail 
both in the dust and radio continuum emission (Omont 
et al. 2001; Beelen et al. 2004) and in the J=5^4, 4^3, 
2^1, and 1— s-O transitions lines of CO (Cox et al. 2002; 
CariUi et al. 2002). With an apparent far-infrared (FIR) 
luminosity of 3 x 10^'^ Lq, it harbors a massive reservoir 
of molecular gas which is the site of active star formation 
with an implied star formation rate of ~ 900 Mq yr~^. The 
CO line emission is resolved into an Einstein Ring with a 
diameter of I'.'S, a direct indication of lensing. The ampli- 
fication factor is estimated to be about 3.5. The data are 
consistent with a disk surrounding the QSO with a radius 
of 2 kpc and a dynamical mass of a few 10^° Mq (CariUi 
et al. 2003). 

The observations of the molecular gas in 
PSS 2322-1-1944 provide one of the best examples to 
date of active star formation in the host galaxy of a 
luminous, high redshift QSO. The strong CO emission in 
PSS 2322-f 1944 makes it a prime target to search for the 
3Pi^3pp transition of [Cl] {vrcst = 492.161 GHz) which, 
at z = 4.12, is shifted into the 3 mm atmospheric window. 
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Fig. 1. Observed spectra of the [Ci]^Pi— >'^Po (this paper) 
and CO(J =4^3) hnes (from Cox et al. 2002) toward the 
z = 4.12 quasar PSS 2322+1944. The black horizontal 
lines show the zero flux density level. The red lines show 
the results of Gaussian + continuum fits. 



Here we report the ^Pi^^Pq [Ci] line in 
PSS 2322+1944. After the detection of [Ci] in the 
Cloverleaf, and, recently, in IRAS FSC 10214 (z = 2.3) 
and SMM J14011+0252 (z = 2.5) reported by WeiB et 
al. (2004), these observations represent the fourth 
clear detection of [Ci] in a high-z source. In this pa- 
per, we assume the concordance A-cosmology with 
Ho = 71 km s-i Mpc-\ = 0.73 and ilm = 0.27 
(Spergel et al. 2003). 

2. Observations and data reduction 

Observations of the ^Pi—^'^Pq transition of [Ci] in 
PSS 2322+1944 were carried out with the IRAM Plateau 
de Bure interferometer (PdBI) in a series of observing ses- 
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Fig. 2. Map of the [Ci] ^Pi^^Pq velocity-integrated 
emission toward PSS 2322+1944. The emission is in- 
tegrated between —400 and +400 kms~^. The white 
cross indicates the optical position at R.A. =23:22:07. 25, 
Dec.=19:44:22.08 (J2000.0). The contours correspond to 
multiples of 3cr = 0.34 mJy/beam. Negative levels are 
shown as dotted contours with the same step. The syn- 
thesized beam of 6'.'8 x 5'.'1 (73°) is shown in the lower left 
corner. 



sions between July 2002 and June 2004. The total inte- 
gration time of the useful data is equivalent to 34 hours 
with 6 antennas. In practice, all the observing time was 
allocated during the summer, when antenna maintenance 
takes place. The total on-sky integration time amounts 
to 116 hours often with 4 or 5 (and sometimes with 6) 
antennas. 

We used the interferometer in the D configuration. The 
combined 3.2 mm data result in a synthesized beam of 
6'.'85 X 5'.'09 at a position angle of 73°. The 3 mm receivers 
were tuned first at the red-shifted frequency of the [Ci] 
^Pi^^Po line, i.e. 96.127 GHz at z = 4.119 and then at 
96.063 GHz for the second half of the data. This is the 
reason why the spectrum is noisier at the very ends of 
the band (Fig. At those frequencies, the typical SSB 
system temperatures were w 150 — 250 K. The water vapor 
ranged between 4 and 10 mm on different sessions. The 
580 MHz instantaneous IF-bandwidth were observed with 
a resolution of 1.5 MHz. 

All data reduction were done using the GILDAS soft- 
wares. Standard calibration methods using close calibra- 
tors were applied. The bandpass calibration was done on 
the quasar 3C454.3. The amplitude and phase calibra- 
tion were performed on 3C454.3 and the nearby quasar 
2230+114. Only data with phase noise better than 40 deg 
were used. The maximum position errors at 3.2 mm intro- 
duced by such phase noise is < 0'.'5. The flux calibration 
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Table 1. Properties of the [Ci] and CO(4^3) lines observed toward PSS 2322+1944. 
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Line 


i^rcst ^obs 


Peak Int. 


Position A«FWHM 


/ 


L' 


L 




[GHz] 


[mJy] 


[kms-i] 


[Jy kms^^] 


[IQi'^Kknis-Vc^] 


[IO^Lq] 


[Cl] ^Pi^^Po 


492.161 96.127 


2.4 


-52±25 319±66 


0.81±0.12 


3.4 ±0.5 


1.1 ±0.2 


CO (4^3)t 


461.041 90.048 


10.3 


-3±10 348±40 


3.83±0.40 


16.3 ± 1.5 


5.1 ±0.5 



NOTE. - The C0(4— >3) data are from Cox et al. (2002). The continuum and the line were fitted separately using an a priori 
line window (between —400 to +400 kms~^) defined from the high signal-to-noise CO J =4^3 spectrum. The luminosities are 
not corrected for the lens amplification (m = 3.5). 



is based on the PdBI primary cahbrator MWC 349. The 
fluxes of 3C454.3 (rcsp. 2230+114) varied from 7 to 4 Jy 
(resp. 4 to 2.4 Jy) during the 3 years observing period. 

A standard calibrated uv table was produced and an- 
alyzed both by direct fits in the uv plane (as this avoids 
the deconvolution step) and by making the deconvolved 
image shown in Fig. |5] Both methods confirm that the 
[Ci] ^Pi^^Pq emission \s i) detected, ii) not resolved by 
our observations and Hi) centered on the optical position 
which is also the phase center of the observations. The two 
latter points allow us to derive the source spectrum from 
the real part of the average of all the complex visibilities. 
The main advantage of this method is to directly obtain 
the spectrum from simple operations instead of applying 
the usual fitting procedures in the uv plane. Any further 
processing (smoothing, fitting) was performed using the 
CLASS software. 

The final [Ci] spectrum at a velocity resolution of 
90 kms~^ is displayed in Fig. ^ together with the spec- 
trum of the J =4^3 transition of CO. The image of 
the [Ci] emission is shown in Fig. [3 The [Ci] ^Pi— >^Po 
fine-structure line is clearly detected at the same red- 
shift as the CO emission with a peak flux density of 
Si, — 2.69 ± 0.31 mJy (including the continuum flux of 
0.31 mJy). Table ^ summarizes the line parameters de- 
rived from continuum + Gaussian fits. 

3. Results and Discussion 

The integrated flux of the [Ci] '^Pi^^Pa fine is 0.81 ± 
0.12 Jykms"^, a factor of 5 lower than the CO(4^3) line. 
A Gaussian fit yields a line width of 319 ±66 km s^^ and a 
line center displaced by — 50kms~^ relative to CO(4^3). 
Due to the weak signal to noise ratio of the [Ci] data, 
the difference in the line position with the higher quality 
C0(4— >3) data should not be over- interpreted. The [Ci] 
line flux implies a line luminosity of Lqj — 3.4 ± 0.5 x 
10i°Kkms-ipc2 or l.l±0.2x 10*^ Lq (see, e.g., Solomon 
et al. 1997 for the definition of the line luminosity). The 
[Ci] luminosity in PSS 2322+1944 is thus about a factor 
2 lower than in the case of the Cloverleaf where L^j = 
6.1 X 10^° Kkms^^ pc^, as derived from the [Ci] velocity- 
integrated flux of Barvainis et al. (1997). 

The continuum emission is detected at 96 GHz with 
0.31 ±0.08 mJy, which is consistent for dust emission with 
the available photometric data at higher frequency. The 
FIR spectral energy distribution of PSS 2322±1944 is wefl 



reproduced with a grey body of temperature 43 ± 6 K and 
a dust emissivity cx jyi-6±o.3 (ggglen et al. 2004). 

The detection of the [Ci] ^Pi—^^Pq emission line 
in PSS 2322+1944 allows us to estimate the mass of 
neutral carbon Mci- Since the ^P2— »'^Pi transition (at 
Vrcst = 809.342 GHz) is not observed, we assume that the 
•^Pi^^Pp transition is optically thin as is the case for the 
Cloverleaf (WeiB et al. 2003). Under these assumptions, 
we can express the mass of neutral carbon as a function 
of L'ci 

Afci = 5.65 X 10-4 I e(23.6/T..) L^j(3p^^3p„) Mq, (1) 

where Tex is the excitation temperature and Q is the Ci 
partition function. Assuming that Tex is equal to the tem- 
perature of the warm dust Tdust = 43 K, the mass of 
neutral carbon in PSS 2322+1944 amounts to Mci = 
4.3 X 10''Mq, or 1.2 x 1O''M0 after correction for lens 
amplification (m = 3.5). Although the excitation temper- 
ature of the neutral carbon can be different than Tdust (as 
observed for the Cloverleaf - see Weifi et al. 2003), the 
derived mass depends only weakly on Tex- For a range of 
[Ci] excitation temperatures from 30 to 100 K, the mass 
of neutral carbon in PSS 2322+1944 would vary from 1.2 
to 1.4 X 10^ Mq. 

Compared to the mass of molecular gas of M^^^ — 
7 X 10^° Mq after correcting for amplification (Cox et al. 
2002; Carilli et al. 2003), the derived mass of [Ci] implies 
a carbon abundance relative to H2 of [Ci]/[H2]« 3 x 10~^, 
indicating near to solar abundances in this high-redshift 
system. Similar values are derived by Weifi et al. (2004). 
This relative carbon abundance is close to the maximum 
value of 2.2 x 10^^^ found for Galactic dense molecular 
clouds with opacities of 4 — 11 mag, a value which does 
not vary within a factor of a few for larger Ay 's (Frerking 
et al. 1989). 

The [Ci] ^Pi — >^ Po and ^P2^^Pi fines are major gas 
coolants. In the Cloverleaf, the ■^P2^'^Pi line is 2.2 times 
stronger than the ^Pi -^^ Pq fine (Weifi et al. 2003). 
To estimate the carbon cooling, we may assume a simi- 
lar ratio for PSS 2322+1944. To compare this to the CO 
cooling, we added the observed CO line luminosities of 
PSS 2322+1944 up to J = 9 - 8, where for the unobserved 
transitions we adopt the prediction of the LVG model of 
Carilfi et al. (2002). In PSS 2322+1944, we find that the 
CO/[Ci] luminosity ratio is 6, as compared to a ratio of 20 
for the Cloverleaf (see Table E))- However, the total mass 
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Table 2. Comparison of the interstellar gas and dust luminosities in Infrared Luminous Galaxies and in the Galactic 
Center 
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C1" 


cott 


C+ 


Lfir 


CO/C 


^ w ~ V / ^ j J-J Ji* 1 JrS. 


C+ /Lftr 


Rof. 






[Lo] 


[Lo] 


[Lq] 


[Lq] 










PSS 2322+1944'"' 


4.12 


1.0 X 10^ 


6.8 X 10** 


< 1.7 X 10^ 


8.6 X 10^^ 


6.6 


9 X 10"'' 


< 2 X 10"" 


[1] 


Cloverleaft") 


2.56 


7.1 X 10'' 


1.5 X 10^ 




4.2 X 10^^ 


21.7 


3.9 X 10"" 




[2] 


Arp220 


0.018 


1.8 X lO'' 


3.9 X 10^ 


1.6 X 10*^ 


1.2 X 10^^ 


2.2 


4.7 X lO"'' 


1.3 X 10"^ 


[3] 


NGC253 


0.0008 


1.6 X lO'^ 


1.6 X 10*^ 


7.8 X 10^ 


1.0 X 10^" 


10.0 


1.8 X 10"" 


7.8 X 10"" 


[4] 


Galactic Center^''' 




6.7 X 10^ 


2.6 X 10^ 


2.6 X 10*^ 


3.9 X 10** 


3.8 


8.4 X 10"" 


6.7 X 10"^ 


[5] 



NOTE. - t Total C l luminosity, tt q,q luminosity up to J = 8. (a) The luminosities are corrected for lensing; m = 3.5 for 
PSS 2322+1944 (Carilli et al. 2003) and m = 11 for the Cloverleaf (Venturini & Solomon (2003). (6) The central 5 x Ideg^ of 
the Galaxy. References: [1] This paper - [2] Weifi et al. (2003); Barvainis et al. (1997); A. Weifi(private communication) - [3] 
Gerin & Phillips (1998) - [4] Bayet et al. (2004) - [5] Fixsen et al. (1999). 



of [Cl] remains somewhat smaller than that of CO (see 
CariUi et al. 2002). 

Compared to the far-IR luminosity of PSS 2322+1944, 
LpiR ^ 8.6 x 10^^ L0 (corrected for lensing), the CO and 
[Ci] cooling represents ~ 10"'' of the far-IR continuum, 
again not very different from the ratio of 4 x lO""* derived 
for the Cloverleaf (Table |2J). The recent search of the red- 
shifted [Cii] fine-structure line in PSS 2322+1944 implies 
an upper limit to the [Cii] line luminosity of 1.7 x 10^ L©, a 
weakness which is typical for high-0 IR luminous galaxies 
(Benford et al. 2004). The CO and C cooling is therefore 
more than half of the cooling due to C+. The [Cii] line 
remains the main cooling line of the gas. 

Finally, the [Ci] luminosity relative to the integrated 
far-IR luminosity is a good measure of the intensity of 
the non-ionizing UV radiation field in galaxies, because in 
photodissociation regions the column density of neutral 
carbon is mostly insensitive to the UV field, whereas the 
far-IR emission is directly proportional to the strength 
of the UV field (see, e.g., Kaufman et al. 1999; Gerin & 
Phillips 2000). For PSS 2322+1944, this ratio is - 3x 10"^ 
indicating a UV radiation field of a few 1000 times larger 
than in the solar vicinity. Both this ratio and the implied 
strength of the UV illumination are comparable to the 
values derived for the other IR luminous galaxies and the 
galactic center listed in Table |21 

As in the case for local starburst galaxies such as 
NGC253, CO is a more important coolant than C in both 
PSS 2322+1944 and the Cloverleaf by about one order of 
magnitude (Table - see also Schilke et al. (1993) and 
Bayet et al. (2004). Similarly, in starburst galaxies or in 
galactic nuclei (including the Milky Way), the molecu- 
lar gas is warm therefore populating the higher CO levels 
which contribute to the cooling. 

The detection of the red-shifted [Ci] ^Pi^^Pq tran- 
sition line in the gravitationally lensed z — 4.12 QSO 
PSS 2322+1944 enables to further constrain the physi- 
cal conditions of the neutral gas and to compare the ma- 
jor line (C, CO and C+) and far-IR luminosities in this 
high-z galaxy. Together with the recent [Ci] detections in 
other sources at high redshift (Weifi et al. 2004), these 
results illustrate the potential of studying neutral carbon 



or species other than CO in high-2 sources, a field which 
will clearly fully develop as soon as more sensitive sub- 
millimeter arrays, such as the Atacama Large Millimeter 
Array (ALMA), will become operational. 

Acknowledgements. We thank the TRAM Plateau de Bure staff 
for their support in the observations. IRAM is supported by 
INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
Estelle Bayet and the referee. Dr. Phil M. Solomon, are kindly 
acknowledged for helpful comments. 

References 

Barvainis, R., Maloney, P., Antonucci, R., & AUoin, D. 1997, 
ApJ, 484, 695 

Bayet, E., Gerin, M., Phillips, T.G., & Contursi, A. 2004, A&A, 
in press 

Beelen, A., Benford, D., Cox, P., et al. 2004, submitted to ApJ 
Benford, D., Cox, P., et al. 2004, submitted to ApJ 
Bertoldi, F., Carilli, C, Cox, P., Fan, X., Strauss, M., Beelen, 

A., Omont, A., & Zylka, R. 2003, A&A, 406, L55 
Bertoldi, F., Cox, P., Neri, R., et al. A&A, 409, L47 
CariUi, C.L., Cox, P., Bertoldi, F. et al. 2002 ApJ 575, 145 
Carilli, C.L., Lewis, G.F., Djorgovski, S.G., Mahabal, A., Cox, 

P., Bertoldi, F. & Omont, A., 2003, Science, 300, 773 
Cox, P., Omont, A., Djorgovski, S., et al. 2002, A&A, 387, 406 
Gerin, M. & Phillips T.G. 1998 ApJ 509, L17 
Gerin, M. & Phillips T.G. 2000 ApJ, 537, 644 
Greve, T., Bertoldi, F., Small, I. et al. 2004, submitted to 

MNRAS 

Fixsen, D.J., Bennett, C.L., & Mather, J.C. 1999, ApJ, 526, 
207 

Frerking, M.A., Keene, J., Blake, G.A., Phillips, T.G. 1989, 
ApJ, 344, 311 

Kaufman, M.J., Wolfire, M.G., Hollenbach, D.J., & Luhman, 

M.L. 1999, ApJ, 527, 795 
Omont, A., Cox, P., Bertoldi, F., McMahon, R.G., CariUi, C.L., 

& Isaak, K. 2001, A&A, 374, 371 
Schilke, P., Carlstrom, J.E., Keene, J., & Phillips, T.G. 1993, 

ApJ, 417, L67 

Spergel, D.N., Verde, L., Peiris, H., et al. 2003, ApJS, 148,175 
Solomon, P.M., Downes, D., Radford, S.J., & Barrett, J.W. 

1997, ApJ, 478, 14 
Solomon, P.M., Vanden Bout, P.A., CariUi, C.L., & Guelin, M. 

2003, Nature, 426, 636 
Venturini, S. & Solomon, P.M. 2003, ApJ, 590, 740 



J. Pety et al.: [Ci] in a quasar at redshift 4.12 5 

Weifi, A., Henkel, C, Dowries, D., & Walter, F. 2003, A&A, 
409, L41 

Weifi, A., Dowries, D., Henkel, C, & Walter, F. 2004, A&A, in 
press 



